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ABSTRACT 


UPDATED POSITION AND ICE VELOCITY 
FOR THE AIDJEX MANNED CAMPS 


A.S. Thorndike and T.O. Manley 


During 1975 and 1976, measurements of sea ice motion were made as 
part of the Arctic Ice Dynamics Joint Experiment (AIDJEX). In an earlier 
report (Thorndike and Cheung, 1977), energy at the inertial period was at- 
tenuated by approximately 50% due to the choice of a low covariance factor, 
Q = 100, during Kalman filtering. 

For oceanographic investigations, the data set proved to be unac- 
ceptable particularly in the summer months when there is a significant amount 
of energy at the inertial period. 

This report is the result of the reprocessing of the original position 
data taken from the four manned camps (Caribou, Blue Fox, Snowbird and 
Big Bear) with a covariance factor Q = 100000. 

The higher Q allows data of the inertial period and longer to pass 
through the filter with less than 1 percent attenuation. 

Because of the large amount of information contained in the data set, 
it was divided into two parts. This report (Volume 1) provides position and 


ice velocity from 11 April 1975 to 17 October 1975. Volume 2 covers the re- 


mainder of the data, 18 October 1975 to 4 May 1976 (Thorndike and Wa 


Accession For 
1980). 
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Introduction 


The objective of the Arctic Ice Dynamics Joint Experiment, as set 
forth in its most general statement, was to reach an understanding of the 
dynamic and thermodynamic interaction between the ice cover and its en- 
vironment (Maykut et al., 1972). One component of the experiment was an 
array of drifting stations at which frequent measurements of the position 
were made (Figure 1). For the present report those measurements have 
been edited using Kalman smoothing techniques to give the best estimates of 
position and velocity. (The raw measurements themselves are not reported 
here, but are available through the Polar Science Center.) The data set 
described here is restricted to the main AIDJEX experiment, which began in 
March 1975 and ended in May 1976. 

This data report represents the most complete documentation of ice 
motion available, but it is limited to one year and to the Beaufort Sea. Par- 
ties interested in earlier data, more widely dispersed in space and time, 
should contact ‘the Polar Science Center, University of Washington, Seattle, 
Washington, which has compiled most of the ice drift data beginning with the 


drift of the FRAM in 1893 to the present. 
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aA full description of the measurement systems used in this program 
is beyond the scope of this report. The list of references includes sources 
of more detailed information (Wasilew and Vivian, 1976; Martin, P., 1978). 
All of the hardware we used was supplied by the Satellite Positioning Cor- 
poration, Houston, Texas. 

Four manned camps -- identified by radio call names Big Bear, Cari- 
bou, Blue Fox and Snowbird -- were tracked using NavSat (Navy Navigation 
Satellite System). The geographical position of each of the four camps was 
measured at unevenly spaced times using Doppler satellite navigation tech-~- 
niques. The raw measurement is the frequency of a signal transmitted from 
the satellite to the ice station. The measured frequency is affected by a 
Doppler shift related to the rate of change in distance between the satellite 


and the ice station, which is itself related to the unknown ice coordinates 


and the known satellite orbit. During each satellite orbit several frequency 


measurements were made and from them an overdetermined solution was found 
for the unknown coordinates. Data recording in the field was on magnetic 
tape with hard copy printout for backup. 

A constellation of six polar orbiting satellites was maintained for this 
system. The orbital period is about 107 minutes. An observer at high lati- 
tudes can track every orbit. Of the 81 potential passes each day, many are 
unusable (because of radio interference between overlapping passes, for in- 
stance). We averaged about 25-30 usable passes per day at each station. 


Single channel (100 MHz) receivers were used to count the Doppler cycles 


for roughly 30-second intervals. The duration of each Doppler counting 


interval was measured independently. The frequency standard used in the 


receivers was stable during a 20-minute pass to about | part in 101}, 


Fix solutions were done in real time at each manned camp. The 
solution used only the Doppler data from the 5 minutes on either side of the 
time of closest approach of the satellite. No additional editing was performed. 
A correction to each Doppler cc 1t was based on the duration of the counting 
interval to account for the random variations in the counting interval. The 
correction also included the fixed timing delay in each receiver. The fix 
solution ignored the effect of the velocity of the ice camps. 

In fact, two receivers with antennas separated by about 100 meters 
were operated at each manned camp. In addition to serving as a backup, 
the measurements at the second antenna were used to define the azimuth of 
the line joining the two antennas. This was possible when both receivers 
produced acceptable fixes from the same satellite pass. Estimates of the 
azimuths obtained at each of the manned camps were reported by Thorndike 


and Cheung, 1977. 
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Preprocessing 
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| The original data tapes from the manned stations measurement sys- 


tems typically contained several "glitches" per tape because of power failures 
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at the camp, or operator error, or some other reason. These tapes were 


transcribed discarding those records which had either a parity error or an 


is nanan ie 


bj incorrect number of words. The tapes which resulted from this process 
were read without error on the CDC 6400 at the University of Washington, 


where the bulk of the processing was done. i. 


Stee 


Three problems had to be dealt with to prepare the sequence of 


measurements at each station for the actual Kalman filtering. First, data 


29) 


recording problems had created occasional gaps in the manned camp data. 
At times the backup teleprinter records from the field could be used to fill 
these gaps. This was done by keypunching the fixes directly from the 
printed records. Second, the data sequence had occasional points which 
appeared to be out of sequence in time. This was sometimes correctable, 
as when the operator had keyed in the wrong date, but usually the situation 
was ambiguous and points out of chronological order were dropped. The 
third and most serious problem involved bad fixes. A glance at a plot of 
the raw time series (Figure 2) illustrates the problem. Because outlying 
fixes have a deleterious effect on the Kalman filtering results, an attempt 
was made to remove them during preprocessing. The algorithm employed 
compares each fix in latitude with the median of the latitude of the 10 fixes 
preceding and the 10 fixes following it. When the differences exceeded a 


preset tolerance of 10 kilometers, the fix was eliminated. The algorithm 
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& was applied twice, with a smaller tolerance of 1 kilometer the second time. 


i The same procedure was then applied for longitude. 


Visual checks were an important part of the procedure. Plots of 


the data points were produced before and after preprocessing. At times it 
seemed from the plots that we were not detecting certain bad fixes, or that 
we were throwing out good ones, and it was necessary to adjust the tcler- j 
ance levels. Edited raw data plots like the one shown in Figure 3 provided 


assurance that the data were ready for filtering. 


For the main experiment, AIDJEX adopted the convention of number- 
ing days consecutively, beginning with day 1 = 1 January 1975 and ending 
with day 500 = 14 May 1976. Throughout this data report, this numbering 


B system is used quite extensively. Appendix 3 provides a conversion table 


from AIDJEX days to calendar days. 
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Kalman Smoothin 


In the smoothing used here an assumption is made regarding the 


motion of the ice. If the state of the ice at time t., is represented by a 


vector X ao containing, say, the position, velocity, and acceleration of a 
point, then the state at some future time is assumed to be partially deter- 


mined by Xn : 


eel Wie (1) 


where $ and © are matricies known from the physics of the situation, and 
Wa +1 Tepresents an unknown random perturbation. In Kalman's formulation, 


which we follow, the random effects are assumed to be Gaussian and White: 


Qifn=m, and 
cov (W_LW_) = (2) 
ete 0ifnéem 


24r-®,  Qualitatively choosing a small value of 


We have chosen Q = 100000 m 
Q is equivalent to assuming that the ice moves in a completely predetermined 
way. (See from equations 1 and 2 that Q = 0 would imply that Snel = 7Xh 

a large value of Q corresponds to assuming that the ice experiences large 
sudden changes in its acceleration.) Thorndike and Cheung (1977) used 

Q = 100 mhr®, The lower value of Q assumed smooth behavior, but was 
unnecessarily conservative. The present value of Q was determined by trial 
and error. <A sequence of raw fixes was used as input to the Kalman filter. 
Several different runs were done using different values of Q. The present 


value of Q appeared to give the best agreement between the output from the 


filter and the raw fixes. 
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The original Q (100 m?nr7®) used by Thorndike and Cheung removed 


approximately 50% of the inertial period amplitude as can be seen in Figure +. 
A constant amplitude sinusoidal function whose period varied linearly with 
time was chosen as the input to the Kalman filter with Q = 100 m?nr~®, In- 
put periods ranged from 0.5 hours (at day 10) to 60.4 hours (at day 15) with 
the inertial period being located at day 10.9. An approximate 50% attenuation 
of the input function at the inertial period can be roughly estimated from the 
graph. 

Figure 5 shows the same input function and the resulting Kalman 
output with Q chosen to be 100000 m?nr®, The inertial period passed through 
the filter with very little attenuation. Another test was then run, this time a 
constant amplitude, constant period (12 hours) sine wave was used to simulate 


the inertial period of the Arctic latitudes. The results indicated that there 


was only a 0.03% attenuation from the original input function. 


For the Kalman reduction, we defined A = ta th and the ; 
matrices 
ease 9 
x = position 2 = Le. i Sade 

velocity 1 A (3) 

acceleration 0 0 0 

r= |1%6 
3 

A7*/2 (4) 
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The dimensions of the matrices are X : 3x1, Q:3x3, 7: 3x1, 9:3x3, 
W:1lxl,andQ:1x1l. 

An identical equation can be written for the y-coordinates. The 
two coordinates are treated independently. 


The measurement Zn is related to X by 


Z = (1 0 0) X,+V 


n n 


; ; ; s 2 
where . a is a measurement error with assumed covariance R = (60 meters)”. 


The response of the measurement and processing system was investi- 
gated as follows. Four hundred sample points were drawn independently 
from a Gaussian distribution with zero mean and variance (50 meters). A 
sampling interval of one hour was assumed. This noise simulated the Nav- 
Sat measurements. The sequence was used as input to the Kalman filter 
procedure. The velocity output from the filter shows two effects. First, 
since velocity is the time derivative of position, we expect to find the spec- 
trum SC) behaving as ws (ud, where S,() is the spectrum of position. 
S,¢~) was assumed to be white; S,() = a* , which implies sy) = Fs 3° . 
The second effect is the tendency of the filter to smooth out high frequen- 
cies. Qualitatively then, we have something like the diagram in Figure 6. 

The actual spectrum of ice velocity is shown in Figure 7. The dia- 
gram illustrates how the measurement errors and finite sampling rate limit 
the accuracy of the velocity estimates. Most of the variance in the velocity 
errors occurs for the frequencies greater than 3 cycles per day. For lower 


frequencies we can estimate the contributions to the variance as 0 - 1.3 cycles 
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Figure 6 - Diagram showing qualitatively, the response of 
the measurement and processing system. 
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per day: g? 33x10 msec”, 1.5- 3.0 cycles per day: g* 3 2x107° 
m2 sec? total variance 0 - 12 cycles per day: 3? z 104m? sec’. 

The general smoothing problem is to estimate X at time th using the 
entire set of observations (2;) i=1....m, and to give the error variance 
of that estimate. The solution equations, due to Kalman, are reproduced 
in Appendix 1. AL positions were converted to an X-Y coordinate system 
prior to actual filtering. This coordinate system along with its governing 
equations are shown in Appendix 2. Typical estimation error variances for 


position and velocity are 40 meters and 1.0 cm/sec respectively. 


15 


Results 


The final data listing consists of five major columns. The first 
column is the time in decimal AIDJEX days corresponding to the observation 
of the position and ice velocity at a particular camp. The remaining columns 
indicate the camps for which data are applicable. Subdivisions under each 
major camp heading are LAT (latitude in decimal degrees), LONG (longitude 
in decimal degrees), and the north and east components of ice velocity in 
cm/sec, NIV and EIV respectively. 

Due to the large amount of information contained in the data set, it 
was divided into two parts. This report (Volume 1) covers AIDJEX days 
101 (11 April 1975) to 290 (17 October 1975). Volume 2 covers AIDJEX days 
291 (18 October 1975) to day 490 (4 May 1976) (Thorndike and Manley, 1980). 

Figures 8- 11 show a more detailed drift track of each of the manned 
camps. An insert shows where the region lies with respect to the north slope 
of Alaska and Banks Island. Asterisks superimposed on the drift track indi- 
cate the positions at integral multiples of 20 days. The beginning and end- 
ing days are noted for each trajectory. A dashed line indicates a long period 
of missing data. The plotted region is 500 x 500 km and is aligned with the x, 
y coordinate system shown in Appendix 2. Figures 8 - 11 were reproduced 


from Thorndike and Cheung (1977). 
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APPENDIX 1 


KALMAN SMOOTHING EQUATIONS 


state vector (n x 1) 
transition matrix (n xn) 


coefficient matrix for random perturbations to physical 
model (nx m) 


white noise driving the physical process (mx 1) 
covariance matrix for W (mx m) 

measurement vector (i x 1) 

matrix relating measurement to state vector (i xn) 
measurement error vector (ix 1) 

measurement error covariance (2 x 2) 


prediction, filtering and smoothing error covariance 
matrix (nxn) 


identity matrix (n xn) 

Kalman gain matrix for filtering (n x i) 
smoothing gain matrix (n x n) 

number of measurements in the data sequence 
index for the data sequence 

time 


The superscript T denotes the transpose overation. The symbol 


< (i j) means the estimate of « at time u given observations (2,, Z,, .... 2). 


a 


The notation here follows Meditch (1969). 


Prediction 


Given observations (Z,, ..., 2.) predict the state at the future 
- 
time tet and find the prediction error covariance. 


Ck +1, ik) = (Kk +1, k) ((k'R), 
P(k +1!k) = 5 (k + 1,k) P(kik) 92 (kK + 1,k) 
+ 7k + 1,k) Qk) f7 (ke = 1,k). 


Filtering 


Given observations (24, Nes Zy) and a new observation ZK +1)? 
estimate the present state at time heed and find the filtering error co- 
variance. 


Lik + if[ k+l) = Vie + 1[k) + K(k + D (Zk +) 
- H(k +1) X(k + 1/k)), 
P(k+1lik+1l) = (I-K(k+ 1 H(k+1)) P(kK+1;k), 


where 


K(k +1) = P(e +k)? (k + 1 (Ak +1) P(e + LIK) AT (k + 1) 


RD) A(5) 


Smoothing 


Given the complete data set (Z,; Beats zy), estimate the state at 
any point the and find the smoothing error covariance. 


K(k kK) +A (kK) (X(K + 1{N) - ¢(K +1; k)) 


= P(k}k+A(kK)(P(k + 1!N) 
- P(k+1|k)) A! (k), where 


P(kik) 3? (k+1,k)P h(k + Lik). 
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Implementation of the Filter Algorithm 

Starting from an initial condition x (0; 0), P(0'0), a prediction is 
made to time the The prediction error covariance matrix P(1!0) is used to 
obtain the Kalman gain K(1). The filtered estimate ((1/ 1) is then found 
by adding to the predicted x (1!0) a quantity which is proportional to the 
difference between the observation Z(1) and H(1)%(1'0). H(1)%(1/0) is 
the value the measurement would have if the true state were actually <(1 0) 
and no measurement errors were made. The algorithm is repeated as each 
new data point is entered until the last filtered quantities { (N | N) and 
P(N |N) are produced. 

The smoothing algorithm, proceeding backward in time, utilizes 
results from the prediction and filtering calculations. The first step is to 
find ( (N-1/N) and P(N-1{N) from equations A(6), A(7) and A(8). No 
direct use of the measurements is made in smoothing. The algorithm simply 
corrects the filtered vector x (kK |k) by a quantity proportional to the differ- 
ence between the prediction vector at time teed and the smoothed value at 
time tee 

In the calculations done for this application, the CDC 6400 computer 
time required to do the prediction, filtering, and smoothing is about 0.5 sec- 
ond per data point, nearly half of which is absorbed in the matrix inversion 
in A(8). Because of the recursive nature of the filter, demands on central 


memory storage space are very modest. 


Reproduced from Thorndike, A.S., 1973, An integrated system for meas- 
uring sea ice motions. In Ocean 73, IEEE International Conference on En- 
gineering in the Ocean Environments, IEEE publ. no. 73 CHO 774-0 OEC, 
pp. 490-499. 


23 


TET OTTERS IN s terre ey 5 


APPENDIX 2 


COORDINATE SYSTEMS 


(km) 
- 1500 ~!000 


Position measurements were made in geographical coordinates 
(latitude north, longitude east). The smoothing operation was 
done in a Cartesian system (x,y), where 


110.949 (90° - latitude) cos (longitude) (km) 
110.949 (90° - latitude) sin (longitude) (km) 
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APPENDIX 3 


CONVERSION TABLE FOR AIDJEX DAYS 
TO CALENDAR DAYS 


For the main experiment, AIDJEX adopted a convention of number- 
ing days consecutively, beginning with day 1 = 01 January, 1975 and ending 


with day 500 = 14 May, 1976. 


In the conversion table, the first column is the AIDJEX day, the 


second is the corresponding day of 1975 or 1976 and the third entry is the 


calendar date. 
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Position and Ice Velocity Data 


The following section comprises the updated positions and ice veloci- 
ties for the four manned camps Caribou, Blue Fox, Snowbird and Big Bear. 


Data span from 11 April 1975 to 17 October 1975. 
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